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Abstract 

The method of integral transforms is reviewed. In the framework of this method 
reaction observables are obtained with the bound-state calculation techniques. 
New developments are reported. 
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I. INTRODUCTION 



The approach discussed in the paper simplifies solving reaction problems. In particular, sub- 
stantial simplifications occur at intermediate energy when many reaction channels are effective. 
No specification of reaction channels at solving the dynamic equations is required, and to obtain 
reaction observables only the bound-state type calculations are to be performed. One avoids 
the consideration of the complicated configuration-space asymptotic behavior of continuum 
wave functions or the complicated momentum-space singularities. In case of perturbation- 
induced inclusive reactions the cross section is obtained in a closed form, while the summation 
of contributions from various final states with the same energy is avoided. The long-range 
Coulomb final state interaction is automatically taken into account in this case. 

In [|T| it was suggested to compare with experiment the Stiltjes transform of response func- 
tions for perturbation-induced reactions, such as the (e, e') reaction. The transform is directly 
obtained with the bound state few-body methods. A qualitative study of the Fourier transform 
of (e, e') response functions was performed at the same time 0|| but with no intention for a 
complete quantitative evaluation. In Q the general procedure presented below was formulated 
for perturbation-induced inclusive and exclusive reactions, and for nonperturbative processes. 
Despite the successful solution of simple models, use of real values of the transform variable 
in that formulation imposed too high requirements to precision of the transform ||. In || 
the real transform variable was changed to the complex one which substantially improved the 
accuracy of the inversion. The application of the method separately to various multipoles also 
increased the accuracy ||. Successful applications to three- and four-nucleon (e, e') and total 
photodisintegration cross sections were given in [0-0] , and the four-nucleon spectral function 
was obtained along these lines [13]. In the response function case, a similar procedure accom- 



panied with the statistical averaging is attempted in condensed-matter physics, see e.g. [[14 
Different from the transforms considered below it is based on the Laplace transform calculable 
with the Monte-Carlo path integration which reliable inversion is a very difficult task. 



II. OUTLINE OF THE METHOD 

Let us consider " E'-overlaps" defined as 

r(E) = J d 7 (X2^ 7 >(* 7 |xi><H£ 7 - E) (1) 

where the state xi (but not necessarily X2) is localized, i.e. its conventional norm is finite. 
The functions \l/ 7 are solutions to the Schrddinger equation: (H — E 1 )'^ 1 = 0, where H is the 
Hamiltonian of a system. The set \l/ 7 is supposed to be complete and orthogonal. The notation 
/ c?7 implies an integration over the continuum states plus summation over the bound states 
in the set. We thus have J g?7|\I/ 7 )(\I/ 7 | = 1. The ^-overlaps (jl]) are the basic quantities of the 
approach, and in this sense they replace continuum wave functions the conventional approach 
deals with. 

The main points are the following. 1. Reaction observables are expressed in terms of the 
_E-overlaps r(E). 2. The latter are computed indirectly[] with the bound state calculation 
techniques. In the next section we explain the first of these points, while in Sec. 4 and 5 the 



1 The functions ^ 7 are not used for obtaining r(E), and continuum wave functions thus do not enter 
the calculation at all. 



second one is explained. Few-body calculations done so far with the method are outlined in 
Sec. 7. In Sec. 6 other non-conventional approaches having some common features with the 
present one are discussed. 

Eq. (HD includes the case of response functions for perturbation-induced reactions. They 
arise if one sets in (1) xi — Oi^o, X2 = O^q, and E = E + e, where are the transition 
operators, and \l/o and Eq are the initial-state wave function and energy: 

R(e) = (l/2)(x + x*), 
x = J d 1 (6 2 ^o\^j}(^y\d 1 ^o}S(E 1 -Eo-e). (2) 

In Eq. (0) \l/ 7 and E y are the final state wave functions and energies, respectively. The 0\ = O2 
case is typical. 

The method proceeds in the following steps ||. First some integral transform with a smooth 
kernel K of the i?-overlap r(E), (or of the response -R(e)), 

J K(a, E)r(E)dE = $(<x) , (3) 

is calculated instead of r(E) itself. To this aim, let us multiply Eq. ([!]) by K(a, E) and 
integrate over E. One obtains 

$(<r) = J d 1 { X 2\^ 1 )K{a,E 1 )^> 1 \xi) 
= Jdj(x2\K{*,H)%){%\xi}. 
Using the closure property one gets 

$(<t) = ( X3 \K{<t,H)xi). (4) 

Eq. (H) may be viewed as a generalized sum rule depending on a continuous parameter. The 
expression @) is evaluated with the bound-state calculation techniques as described in Sec. 
4. At the second step, Eq. (|3|) is considered as an integral equation to invert and thus one 
gets r(E) from $(cr). In case (0) this provides a final result since the response function is an 
observable quantity. In the general case reaction observables are expressed as some quadratures 
in terms of i?-overlaps (P thus obtained, see the next section. (The reader interested only in 
the calculation of response functions may skip the next section except, perhaps, Eqs. (|TT])-(|T(J). 
Note also that in the response function case the transform variable is changed to e = E — Eq.) 
In the non-relativistic case, all the considered states, or operators may be viewed as "intrinsic" 
ones. 



III. REACTION OBSERVABLES IN TERMS OF £ OVERLAPS 

First let us consider the case of transitions to continuum induced by a perturbation O. The 
matrix element 

M fi = (^J{E)\d^ Q ), (5) 
is to be calculated where \l/o is a localized initial state wave function, while tyj belongs to 



continuum. We set 15 



tyj(E) = (j)f(E) + (E — H — ir]Y l (j)f(E) . (6) 



Here 0/ is either the "channel" function that consists of a product of wave functions of frag- 
ments pertaining to a given reaction channel^, or the sum over such channel functions corre- 
sponding to (anti)symmetrization [|15|]. In the normal non-relativistic case one may use intrinsic 
channel functions that are angular momentum coupled products of intrinsic wave functions of 
fragments and of the relative motion function with a given orbital momentum. The functions 
4>f are 

fo=(H- E)<t> s . (7) 

Using Eq. (|5|), we rewrite (|5|) as 

M fi = (<f> f (E)\6%) + (ME)\(E -H + i-ny^Vi) . (8) 

The first, Born, term in (|8]) is computed directly. The second term takes into account the final 
state interaction. We present this term in the form 

Y,(E - E n )- l F fi (E n ) + (E-E' + i V y l F fi (E')dE' 

n JE th 

roo 

= J2(E - E n )- x F fi (E n ) - mFfi(E) + P (E - E')- l F fi {E')dE' , (9) 

n JEth 

where the form factor Ffi is defined as 

Efi(E') = J d 1 (4> f {E)^ 1 )^ 1 \6^ i )8{E 1 - E') . (10) 

Eq. ( |IUD is a special case of Eq. (|I|), and the reaction amplitude (|5p is thus expressed in terms 
of a simple quadrature @ over a i?-overlap ||. 

The above formulation is applicable when the long-range Coulomb interaction of fragments 
in a final state is neglected. In order to take into account this interaction one can proceed 
as follows. To simplify notation we omit (anti)symmetrization of the channel wave functions. 
(Changes to include (anti)symmetrization are obvious.) Let us write H = Tf + V} where Vf is 
the interaction of fragments in the final state. Let us subtract from Vf a part Uf that includes 
the long-range Coulomb interaction: Vf = Vf + U9 ' . Let us denote <pj the "plain wave plus 

ingoing wave" continuum wave functions of the Coulomb-like Hamiltonian Hf = Tf + U^. 
Instead of Eq. (§ let us use (cf. |S|, Chap. 5, Eq. (108)) 

tyf~ = tpj + (E — H + irj^VfipJ. 

One can proceed as above with the replacement of the channel function 0/ by the "Coulomb" 
channel function ipj, and 0/ by VfipJ, or, more generally, by (H — E)ipJ if <pj is 
(anti) symmetrized. This procedure requires constructing ipj which can be done exactly when 
the Coulomb interaction in only one pair of fragments is present. 

The arguments similar to those presented above are valid if one replaces O^o with any 
other localized state. Consider, for example, a breakup reaction amplitude in the framework 
of the Glauber approximation. It is given by the matrix element of the (approximate) fixed- 
center scattering amplitude between the initial and final state of the target. We thus encounter 
again a matrix element of the form {^>f\x), ^ f and \ being the final state wave function and 
a known localized function, respectively. Proceeding as above in connection with Eq. (|5p, and 



2 The quantity (0) equals Vf(j)f in this case. 



putting x instead of O^/q one replaces the calculation of the continuum wave function tyf by 
that of a .E-overlap. 

In Ref. |l3j the spectral function of an A-body system is expressed in terms of an (A — l)- 
body .E-overlap. Below the formulas are presented in more detail. The spectral function 
P(k, E) represents the joint probability of finding a nucleon with momentum k and a residual 
subsystem of A — 1 nucleons with energy E in the ground state at rest. The energy E is 
calculated with respect to the ground state energy: 

Jo 



P(k,E) = (2J + l)- 1 EE/ dfdK o \(i:f K l,k,s z t z \^(J M ))\ 2 5(E ~ (Ef- 1 ~ E^)). 

(11) 



M =-J s z =±l/2' 



A-1 



Here s z and t z are the third components of the nucleon spin and isospin, Ej = Mj — 
(A - 1)M, Eq = Mq - AM, where Nlf~ l and M A are the masses of the final residual state 
and the initial state, and M is the nucleon mass. The functions V^Kq f° rm the complete set 
of eigenstates of the A — 1 subsystem with the total momentum Ko: 

V'/Ko = / dki . . . dk A ^j-\k u k A _ x )5(E k, - K )|k! . . . k A _ x >. (12) 

J i=i 

The ground state can be written as 

A 

/(P)dP J rfk 1 ... c /k A ^(k 1 ,...,k A )5(Ek i -P)|k 1 ...k A ), (13) 

i=i 

where /(P) represents a narrow peak concentrated in the vicinity of P = 0, and J |/(P)| 2 dP = 
1. In Eqs. (|l2j), (0) and below we omit spin and isospin variables. 

In terms of the wave functions entering (|I2"D, ( jl3|) the spectral function can be written as 



P{k ) E) = {2J, + iy 1 E fdfSiE-iEf-'-E^Wfik), (14) 

M =- Jo 



W f (k) 



,. A-1 

\dk x ... dk A ^' l *{k x , k A _ x )^(ki, . . . , k A _i, k A )5(E + k) 

i=l 



(15) 



with k A = k. 

For the A — 1 subsystem and for the A system in Eq. (15]) let us perform orthogonal 
transformations to Jacobi momenta 



7ri, . . . ,7r A _ 2 ,n A -i = (A - 1) 1/2 E k i! and 7ri,...,7r A _ 2 ,7TA-i,n A = A 1/2 E k 

1=1 1=1 

respectively. Here 7r A _ x = [(A - l)/A]^ 2 \k A - (A - l)" 1 E^ 1 kj. Eq. (|TJ) takes the form 
W>(k) = [A/ (A - 1)] 3/2 | dm . . . dnA^ipf- 1 *^, txa-2] -k)^(7n, • • • , vr A _ 2 , 7r A _ i; 0) 



(16) 



with 7r A _! = [A I {A - l)] 1 / 2 k. Here <pj~ l = (A - l)" 3 / 4 ^/ -1 , = ^~ 3/ Vo > and 
J dm... d-KA-iCpjr^tpf- 1 = 6{f - /'), /dTri . . . dTTA-il^r = 1. 



The last arguments of the functions (p^ 1 and ip^ are total momenta of the A — 1 subsystem 
and A system, respectively. In the non-relativistic limit y5/ _1 and tp^ cease to depend on 
these arguments. From Eqs. fll4]) , (|16|) one sees that the spectral function is a parameter k 
depending P-overlap in the subspace of A — 1 particle. 

The amplitudes of nonperturbative reactions have been expressed in terms of P-overlaps 
in Ref. Q. Usual processes with two fragments in an initial state and any number of fragments 
in final states were considered. Similar to Eq. (H) we shall use the functions 

4>i — (H — E)<f>i <f)f = (H — E)(j)f , (17) 



where <pi and <pf are defined in the same way as <pf in (|7|). We rewrite the expression |15| for 
the i — > / transition T matrix in terms of 4>ij'- 



Tfi(E) = Tfr + (MEW -H + iv)-%(E)) , (18) 



where 



Tfr = (hrnuE)) = (me)\me)) 



is the Born term. 

We need to calculate the second term in Eq. (|T8| ) representing the final-state interaction 



correction. Similar to Eq. ( [TOD we define a form factor 

Efi{E, E') = J d-y(ME)\%)(%\h{E))6{ET-E'), (19) 



and, finally, as in Eq. (|9|), we express the contribution in Eq. ( |T8"| ) we consider in terms of 
the quantity (0) . Again, Eq. fllTf ) is a special case of Eq. ([[]). Because there are two 
colliding fragments in the entrance channel the function 0, is localized. The method is directly 
applicable when the long-range Coulomb interaction is neglected. 



IV. EVALUATION OF THE TRANSFORMS AND THE CHOICE OF THE 

TRANSFORM KERNEL 

Having expressed reaction observables in terms of E'-overlaps, now we consider the calcu- 
lation of the latter quantities using Eqs. (|3|), ([|) of Sec. 2. One possibility is the following. 
Consider some complete set of localized states and denote P/v the projector onto the subspace 
of the first iV states in the set. As iV tends to infinity at some conditions one has 

(X2\PnK((t,H)P n \xi) - ( X 2\K(a,H)\ X i) (20) 

uniformly with respect to a. Therefore Eq. (d) can be approximated with the left-hand side of 
Eq. (|2"0"|). Correspondingly, if one denotes by -ip^ the states that diagonalize the Hamiltonian 
in the subspace determined by the projector P/v and by E^ the corresponding eigenenergies, 
then one has 

N 

= EM^W^XiW'iEW). (21) 

71=1 

This expression can be used as an approximate right-hand side in the integral equation (^). 
Thus, the input to the integral equation is obtained with the bound-state calculation technique. 



Eq. ( pO|) is valid when both xi an d X2 are localized, and the kernel K is bounded uni- 
formly with respect to o: \K(a,E)\ < C. If the latter is valid then the state K{a,H)x\ is 
localized along with xi- Indeed, the norms of these two states are / dj\K(a, £' 7 )| 2 |(^/ 7 |0i)| 2 
and / (i7|(\l/ 7 |0i)| 2 , respectively. Since the latter is finite, the former is finite as well. Taking 
this into account, one sees that the quantities {4>2\QnK4>\), (4>2\KQn4>i) 1 and (4>2\QnKQn4>i) 
tend to zero at N — > 00, Q N being 1 — P^. This gives Eq. fl2~ip . 

The exact solution with the right-hand side (^Tj) is 

jv 

r N {E) = Y,(xM N) )(ri N) \Xi) 5(EM - E) . (22) 

n=l 

While $jv(cr) — > $(cr) as A/" tends to infinity, the amplitudes entering Eq. (|22D are rather 
chaotic in the continuous spectrum region i£ > i£ th , and r N (E) thus does not tend to the 
true r(E) in the normal sense. However, if approximate solutions to the integral equation are 
sought for in the class of smooth functions then the sequence of the approximate solutions 
r smooth^^ w -jj ^end to the true r(E) as N goes to infinity. (See also the next section.) 

In other words, if one averages both r(E) and tn(E) over E using a smooth kernel K(a, E) 
then the corresponding " averaged ^-overlaps" $ and $ at should be close to each other provided 
that N is high enough. If K effectively vanishes for E values beyond some range AE (depending 
on cr), the case most interesting for us, one can expect that the two averaged ^-overlaps will 
be close to each other when a sufficient number of levels E^ fall into AE. Then one can use 
the averaged discretized i?-overlap $at(o") to reconstruct the true i?-overlap r(E) via solving 
Eq. (H) in an appropriate way. 

The procedure is applicable with any smooth K(a, E), say K = exp[— (a — E) 2 /a^}. On the 
other hand, it is restricted to the case when both xi an d Xi are localized. Besides, the complete 
diagonalization of a Hamiltonian is not manageable for larger systems. Let us consider the 
evaluation of the transform with the help of Eq. (|4]) without the above restriction and without 
having recourse to the diagonalization of a Hamiltonian. Unlike the considerations above, this 
evaluation will be done for the kernels of a special form. We set first 

K(a,E) = (E + a)- 1 . (23) 

According to Eq. (||), this leads to the transform 

$( C r) = ( X2 |(F + a)- 1 Xl ). (24) 

Eq. Q2"3D can be written as 

$(a) = (x 2 |*(cr)) , (25) 
where ^ is the solution to the Schrodinger-like equation with a source, 

(H + a)V = X i. (26) 

The a values we use are such that — a lie apart from the spectrum of the Hamiltonian so that 
K is bounded. Due to this reason \1/ is localized along with xi as h is explained above in 
connection with Eq. (^0|). Due to the same reason there exists only one localized solution to 
Eq. (PB|): the difference of any two localized solutions would satisfy the homogeneous equation 
that has no solutions yet. 
One can also use 



K(a,E) = \(E + a*)(E + a)]- 1 



(27) 



We consider here the case of the response function (|2]) and replace E with e = E — Eq in Eq. 
(^7f). According to Eq. @ the transform is 

$((t) = (1/2) [(6 2 *o|(# - + ct*)- 1 ^ - £ + a)- 1 ^!^) + cc] . (28) 

The quantity (|28|) can be calculated as 

= (1/2) [<* 2 (o-)|*i(o-)> + (y 1 (a)\y 2 (a)}] , (29) 

where ^ are solutions to 

(H-E + <T)V 1 = diV Q (H-E q + <t)* 2 = O 2 ^ - (30) 



The first of Eqs. (|30|) is, of course, Eq. ( pq) in the different notation. From Eqs. ( p9|) 
at Oi = O2 it could be inferred once more that are localized. The general method we 
discuss has been given in Ref. [|J with use of the kernels (^3]) and ( p7|) for real a values. The 
corresponding integral transforms are called the Stieltjes and generalized Stieltjes transform 
|16|| . We also note that if Eqs. (^) are solved with the help of an expansion of \I/ over a set of 



known functions then e.g. the quantity fl28|) thus obtained coincides with that in Eq. ( pT|) with 
the kernel (p7|) provided that the diagonalization of a Hamiltonian is done in the subspace of 
retained basis functions. 

So far we chose the kernel in a way that the evaluation of the transform, Eq. (|J), is possible. 
Now let us discuss the choice of the kernel from the point of view of the stable reconstructing 
the response from the integral equation (0). At the same accuracy in $, the inversion of the 
transform (^) will be the most stable if the kernel is chosen as " narrow" as possible. Indeed, 
at integrating in Eq. (|^), variations of r occurring at intervals AE that are smaller than the 
range of the kernel may smear out and even tend to cancel at some a. It then may be hard 
to reconstruct such changes in presence of numerical inaccuracies in $. The kernels ( p3|) and 
( p7D with real a do not possess a "finite" range and thus are unfavorable from this point of 
view. Use of complex o = —07? + 207 in Eq. (|27|) has been suggested in the response function 
case H to increase the stability. The kernel then takes the "Lorentz" form, and 



E tt: %~^ + u' ff. 3 * = »(*)■ ( 31 ) 



Rn + f°° R(e) 

{o-r - e n ) 2 + erf Je th (a R - e) 2 + aj 

Here the contribution from discrete levels is written down explicitly. For a stable reconstruction 
of the continuous part of the response it is advantageous to use mostly cr^ > e t h, and 07 should 
be chosen "sufficiently small". In some sense 07 plays the role of the energy resolution: any 
two responses that noticeably differ from each other at Ae intervals not much smaller than 07 
lead to noticeably different transforms $ and thus can be discriminated even in presence of 
numerical uncertainties in $. But in order to distinguish between responses differing from each 
other in regions small compared to 07 a higher numerical accuracy in $ than that required in 
R might be necessary. 

The procedure similar to that used above to calculate the Lorentz transform in the response 
function case can also be used for other Z?-overlaps provided that in Eq. (JJ) both xi an d Xi 
are localized. The following procedure can be suggested in the general case. It is applicable 
for breakup reactions with three or more fragments in the final state when only Xi is localized. 
Ons writes down the kernel (^7]) as 

K(a, E) = (2ia I )- 1 [(E + a*)' 1 ~ (E + a)" 1 ] , (32) 



and, similar to Eq. (p5|) , one has 



$(a) = (2^ / )- 1 [(x 2 |^ , )-(x2|^)], 
where and \& are, respectively, the localized solutions to the equations 

(H + a*)& = Xi (H + <r)* = Xi- (33) 

In the response function case one may replace E by e = E — Eq, and H by H — E in Eqs. 
(0) and (H), respectively.0 

The kernel (^) with complex cr = — o R + 207: 

tf(<7, E) = [{E - a R ) - iaj]/[{E - a R f + a]} (34) 

is also presumably good from the point of view of the inversion stability. 

We note that no information on exit channels of a reaction is required in order to solve 
the above listed dynamic equations. In the response function case, the long-range Coulomb 
interaction manifests itself only in Eqs. ( |30D and does not cause any problems. 

To calculate the response function, Eq. (Q), it is convenient to use the multipole expansion 
of the transform. This allows taking into account contributions from various projections of 
angular momenta in a general form. We perform the calculation in the reference system where 
the initial state \l/o is & t rest and carries a given total spin J and its projection M . It will be 
denoted ^j m here. We set 0\ = 2 = O in Eq. (0) and consider the response R averaged 
over Mo along with the corresponding transform $: 

R=(2J + l)- 1 £ R(M ) $ = (2J + 1)- 1 £ $(M ). (35) 

M =-Jo Mo -Jo 

We expand the operator O in a sum of irreducible tensor operators Oj m : 

O = ^^djmOjm, (36) 

jm 

where the expansion coefficients. We define the partial right-hand sides 

QjM = ^2 ^jmJ M Ojm^J M , 
m+M =M 



(H-E + a)¥ JM = qUi (37) 



the partial solutions ^j M : 



and the partial transforms: $jj = (q J JM \^jm), or $jj = (^jm\^jm) f° r the case of Eq. 
and Eq. fl2~7|) , respectively. The latter quantities do not depend on M. (These quantities are 
the corresponding transforms of the components Rjj of response functions which correspond to 
transitions to the continuum states with given JM induced by the Oj m multipole operators.) 
One has 

O^J M = ^2 a jrnCjmj oMo qj M ^ = a jm,Cj m j oMo ^ j M . 

jmJM jmJM 



3 It is easy to conclude from Eqs. (|32|), (|33|) that in the response function case the "squared Lorentz" 
transform K 2 (a, E) = [(E — Eq + o~*)(E — Eq + <t)]~ 2 is also calculable in a simple form. The kernel 
K 2 is more sharp peaked than K which facilitates the inversion. 



Putting these expressions into Eq. fl25|) or (^), substituting this into the second of Eqs. (p5|), 
and then performing the summation of the products of Clebsh-Gordan coefficients one finally 
gets 

$ = (2J + l)- 1 £(2 J + l)^j J2(2j + l)- l \a jm \ 2 . (38) 

jj m 



We also note that in Eq. ( 

{VjmWjm) = (2J + l)- 1/2 (^||Oj||^ Jo ) 
where <£>jm is a (basis) state with given J, M values, and the right-hand side includes a reduced 



matrix element (defined as in ||17|| ). 

If O(q) is the charge operator then one has in Eq. (|36|) dj m = Y* m (q), and the sum over 



m in Eq. fl38|) equals (Att) l . Let us consider the transversal current operator J*(q). Eq. ([36] 
should be modified to 



J t (q) = £Y**(q)0> 



^A* i-^\AX 
jm 
jmX 



where Y^ m are vector spherical harmonics ||17|| , and the A superscript signifies "electric" or 



"magnetic". The quantities defined above acquire the additional superscript A, and Eq. ( p8|) 
takes the form 

jJX 

Here the orthogonality property of vector spherical harmonics was used. 

It was implied in the formulas of this section that the transition operators 0\$ in Eq. (0) 
do not depend on the energy transfer e. If a dependence on e occurs it leads to a dependence of 
Oi 5 2 on H after performing the closure. Then the above methods for calculating the transform 
should be modified. In most of the cases the e dependence emerges only as a dependence of 
nucleon form factors on e. The nucleon form factors can be factored out from a calculation, 
and the formulas of the above type then become applicable. Often it is sufficient to divide out 
only one general factor. For single-particle nuclear currents this takes place if one neglects the 
e dependence of the form factor ratio G^/G P E or uses the similarity of the form factors G\ and 
G P M . In case of zero isospin nuclei there is no interference terms between isovector and isoscalar 
components of the transition operator, and the same can be achieved without approximations 
by dividing out the isovector and isoscalar nucleon form factors. 



V. INVERSION OF THE TRANSFORM 

One could compare the calculated transform of a response function with the transform of 
experimental data on the response function |[T|.[TB|.[T^] , and this does not require inversion of 
the transform. However at such a comparison physics might be obscured by the fact that 
parts of different nature of a response are mixed up in the transform. In addition, with 
a given approximate Hamiltonian, often the comparison of a low-energy part of a response 
with experiment is only sensible. At the transform level, the comparison is then obscured by 
high-energy contributions. (These difficulties are less pronounced in case of sharply peaked 
kernels. But for such kernels inversion is also an easier task.) Therefore performing inversion 
is preferable, and below we consider procedures for solving Eq. <^). 



First we disregard discrete contributions to the solution and thus consider the integral 
equation 



/ K(a,E)r(E)dE = <$>(a). (39) 



Afterwards we shall remove this restriction. In the following we assume that Eq. (BSfl has only 
one solution. This is obviously true for the kernel (|23|). Let us demonstrate that this is also 
the case for the Lorentz kernel. The difference Ar(E) of two possible solutions would satisfy 
the homogenous equation of the form 



dx[(x - vY + /3 2 ]- L Ar(x) = 



This leads to 



roo roc 

/ dx dy[{x-y) 2 + (3 2 }- lAr ( x ) Ar (y) = °- ( 40 ) 

J a J a 

Using the representation 

{x 2 + (3 2 Y l = (2/3)" 1 / dkexp(-[3\k\ + ikx), 



(3 > 0, one can rewrite the left-hand side of Eq. (|40| ) as 

(2/5)" 1 / dkexp(-P\k\Mk)\ 2 (41) 



with 

tp{k) = I dx exp(ikx)9(x — a)Ar(x) 



The quantity fl41|) can equal to zero only if (f(k) = that leads to Ar = 0. 

In the calculations performed so far [[4|-|i"3|| the inversion procedure was the following. The 
solution was sought for in the form 

N 

r(E) = Y,CnME,a) (42) 

n=l 

where Xn are known functions forming a complete set. They include non-linear parame- 
ters a. If one substitutes the expansion (^) into the left-hand side of Eq. (^) one obtains 
J2n=i c n L Pn(cr, a) where (p n are transforms of the basis functions. The expansion coefficients c n 
and the parameters a are obtained from the best fit requirement 



A" 

E 

k=l 



N 



n=l 



min. (43) 



At fixed a this leads to a system of linear equations for the expansion coefficients. (At large 
iV this system may become ill-defined. Instead of solving it, the SVD algorithm for a direct 
minimization in ([4"3"D, see e.g. ]2Uj, may be applied.) 



In the Lorentz case, see P-|l"3H the values a k = + i<Jj in ( f43D were chosen such that 
the o"^. points covered the interval < e < e max where e max is such that R(e > e max ) is 
already very small. The <j\ points covered also some interval with e m j n < e < e#, aiming a 
better description of the low-energy part of the response. The 07 values were chosen taking 



into account the following. When 07 is tending to zero K(a,e) — > (7r/|<7f |)<5(<7r — e), and 
$(cr) — > (n/\ai\)R(aji). Hence for sufficiently small 07 the required relative accuracy in R is 
the same as in <3>, and no additional uncertainty in R arises due to the inversion procedure. On 
the other hand, the smaller 07 is the harder is to achieve this accuracy at solving the dynamic 
equation ( |30"D with the bound-state type methods. Indeed, at 07 — > the scattering regime 
is recovered. For few-nucleon systems use of 07 values which are comparable to the widths of 
response functions and fall into the range between several MeV and 20 MeV was found to be 
the most convenient.^ (Of course, the final results should not depend on 07 that can serve as 
a possible check.) Values of a such that a = — ctr + £07(07*) with 07 increasing as cxr increases 
were found to yield good results, in particular. 

Solutions to general equations of the form (|39|) are known to be unstable with respect to 
high frequency oscillations. Indeed, if r(E) is the solution with an exact $ then, for example, 
r a (E) = r(E) + X sin Et is the solution with the right-hand side 

$ a (cr) = $(cr) + A / dEK(a, E) sin Et. 

JE th 

At any A the variation $ a — $ of the right-hand side becomes indefinitely small at sufficiently 
large t. However, the corresponding variation r a — r = A sin Et of the solution may take values 
as large as A. 

Because of the above instability one should not seek for the exact solution of the equa- 
tion with a given approximate right-hand side. Instead, as it is well-known, a regularization 
procedure should be applied at finding the solution. The regularization suppresses a quickly 
oscillating component in the solution with an approximate right-hand side but do not influence 
much a slowly varying component in the solution. If r a and r are solutions with right-hand 
sides $ a and $, respectively, and the difference $ a — $ is small then only a quickly oscillating 
component in the difference r a — r may have large amplitudes. Due to this the regularization 
guarantees the closeness of a regularized "solution" to the exact one. 

When one looks for a solution in the form (^2]) the number N of retained functions plays the 
role of a regularization parameter. If a right-hand side $ a is accurate enough then, normally, 
the accuracy in the solution r(E) first increases when N increases, and the results look as 
being "convergent" with respect to N. But at a further increase in N the "solution" obtained 
inevitably acquires unphysical oscillations. The higher accuracy in $ is, the higher N values 
are at which the oscillating regime starts, and the higher accuracy in r achieved at lower N 
values is. 

Stability of inversion results with respect to N served both for choosing N values and as 
the main quality criterion of an approximate solution. To check the results, responses obtained 
with different sets of basis functions were compared to each other as well. (Such a comparison 
implies that the stability ranges AN exist for each of the sets). Stability with respect to an 
accuracy in the input $, was, of course, also checked. The quality criterion related to sum 
rules was applied in addition. The following sum rules should approximately be satisfied by 
the solution r(E): 

/*oo roc A 

/ r{E)dE={<h\<h), / Er{E)dE — (0 2 |#|0i) ■ (44) 

JEth JE th 

Incorporation of the known threshold behavior of r(E) into basis functions entering Eq. 
( fj2|) increased the accuracy of inversion. Other features of the solution such as positions, or 



For e.g. three-cluster Li response convenient 07 values proved also be comparable to the width 
of the peak. 



widths, of narrow resonances (if known) should similarly be taken into account for the same 
purpose. 

In some cases the accuracy can also be enhanced in the following way M. Let us consider, 
for example, the longitudinal response function (Q) for the (e, e') reaction. It may behave in 
a rather peculiar way at small e. The origin is that the lowest I multipole contributions Ri to 
R — J2iRi have maxima in the threshold region while other multipoles exhibit maxima in the 
quasielastic peak region. It is helpful to apply the whole procedure separately to those lower 
Ri and to the sum of all other Ri. The reason for an increase in accuracy is that with use of 
simple functions ip n in Eq. (E2f) one can better describe the behavior of those pieces of the 



response function at sufficiently low N values than that of the total response function. 

The accuracy of the inversion results has been assessed in Ref. ||) using the two-nucleon 
d(e, e') longitudinal response as an example. An approximate Lorentz transform with a several 
per cent oscillating error was taken as an input. In Fig. 1 the response obtained via inversion 
of the transform at 07 = 5 MeV is shown (dashed curve) along with the exact response (solid 
curve) calculated in the conventional way. The quality of the inversion is excellent. 

Regularization procedures different from that discussed above to solve the general integral 



equation of the first kind were applied in the literature, see e.g. p3|-|22|j. In our case the 



following one, see pT|J2^] , seems to be promising. Let us use the quantity 

1/2 



p(/i,/ 2 ) = / \h{a) - f 2 {a)\Ua 

as a measure of proximity between fi and f2- Let $ a (o") be an approximate right-hand side 
used as the inversion input. We suppose that 

p($ a , $) ~ 5. (45) 

The procedure requires estimating 5. For this purpose one can proceed as follows in our case. 



Let v be some parameter (e.g. the number of basis functions at solving Eqs. (|26p, (|30|)) that 
determines an accuracy in $ a at a given a value: 

= /(!/), $(<r) = /(oo). 

One can estimate the difference j[y) — /(oo) and hence the quantity (^) performing the 
calculation at various v values, then approximating the results with some analytic formula, 
such as e.g. j{y) = /(oo) + CV~ 7 , and fitting the parameters of the formula that include 
/(oo). It is supposed that inversion results are not sensitive to a precision in 5 which normally 
is the case. 

There is no sense to solve Eq. (|5D|) with an accuracy higher than 5. Therefore one seeks 



for the solution in the class of functions r(E) satisfying the condition 

p(Kr, $ a ) < 5 (46) 



where Kr is the short notation for the left-hand side of Eq. ( j39|) . We consider <3> a such that 
X>wT \^a{c)\ 2 da > 5 2 , so that r = does not belong to the class (|46D of functions. One wants 
to find an approximation to the true r(E) that is stable with respect to small variations in $. 
To this aim, the functional 

Q[r) = r \q(E) \r(E)\ 2 +p(E) \dr/dE\ 2 ] dE (47) 



is used, where q and p are non-negative functions. The function r(E) from the class ( f46"D 
at which the functional Q[r] takes its minimal value is chosen as the approximate solution. 



(The minimal value does exist.) If one adds a quickly oscillating component to the true r(E) 
then the contributions to Q[r] from both terms in the right-hand side increase. Therefore 
the function r(E) giving a minimum to Q[r] belongs to the "smoothest ones" among those 
satisfying Eq. fl46|) . (In Eq. fl4"7|) the simplest choices q = l,p = 0, or q = 0,p = 1 are, in 
general, reasonable.) It can be shown |22| that if one takes a sequence of 5 values tending to 



zero then the corresponding sequence of approximate solutions found in the above way would 
tend to the exact r(E). Thus at sufficiently small 5 closeness between the approximate r(E) 
constructed as described here and the exact r(E) is guaranteed. 

It is easy to show [22| that r(E) giving a minimum to Q[r] satisfies the condition pB| ) with 



the equality sign []: 

p(Kr, $ a ) = 5. (48) 

Therefore one can seek for the minimum of Q[r] at the additional condition (fJHJ). This is a 
classical variational problem which can be solved using the Lagrange multiplier method. One 
forms the functional 

p 2 (Kr,$ a )+fin{r], (49) 



finds its minima r M at given /i values and, finally, fixes fx from the equation ( f48|) with r = r M . 
One can show [f22] that the corresponding solution fx = fx(5) to this equation necessarily exists. 

In our case it is convenient to seek for minima of the functional ( |I9|) using the expansion 
of Eq. ( |3^ ) type. At given fx this leads to a system of linear equations for the expansion 
coefficients. If one puts fx equal to zero this system of equations turns to that following from 
Eq. (fj3|). In contrast to the case of Eq. (|43|), in the present case one can retain as many basis 
functions in the expansion as desired. As above to increase the accuracy of the inversion one 
should incorporate known features of a solution into basis functions. 

Since the E'-overlap ([!]) may include discrete contributions ~ 5(E — E n ), the sum of the 
form J2 n K(ct, E n )r(E n ) is present in general in the left-hand side of Eq. Q39T). A good strategy 
is to calculate this sum beforehand which leads to Eq. (|39|) with a modified right-hand side. 

It is convenient to calculate the required r(E n ) from the equation of (|30|) type 

(H - E)>$ = xi 

with real E ranging in the vicinity of the level E n . Then at E — >• E n 

{ X 2\4>(E))^(E n -Ey l r(E n ). 

Another possibility applicable in the framework of the method ( fl~3|) is to solve the corre- 
sponding equation of Eq. ( |3lD type directly with inclusion of the basis functions S(E — E n ) 
into the expansion (f4~2|). 



5 This can be seen as follows. Suppose that r = ro gives a minimum to fi[r], and p(Kro, $o) < 6. 
Then there exist surroundings (p_e(^i^o) < e) of ro such that p(Kr, 3> a ) < 5 for any r belonging 
to the surroundings. For the function r = jtq, \j\ < 1, belonging to these surroundings one has 
f2[r] = 7 2 il[ro] < fi[ro]. Therefore the conditions fi[ro] = min and p(KrQ, $ a ) < S are incompatible. 



VI. SOME OTHER NON CONVENTIONAL APPROACHES 



In Refs. p3|j24|1 the following method for calculating elastic scattering and ionization pro- 
cesses in atomic physics was suggested and used. One calculates the T-matrix ( |18|) for a set 
of complex energies E and then performs an extrapolation to real E values. This leads to 
the dynamic equation of Eq. fl26|) form at a sequence of 07 values tending to zero. Such a 
method probably requires more labour to attain accurate results then the method of integral 
transforms. In contrast to what was said above concerning Eq. (f42|) , within such a framework 
one cannot use the known information about a solution in order to improve accuracy. Probably 



this is the reason for unsatisfactory results obtained in |24| in the threshold region. Besides, 
while one needs the extrapolation 07 — > 0, difficulties in obtaining accurate solutions to Eq. 
(f26|) with bound-state methods increase as 07 decreases. 

In Ref. [25] the following method for calculating response functions is suggested. One can 
represent the 5 function in Eq. (|2|) as — (7r) _1 Im(.E — E 1 + ir])' 1 , and, correspondingly, the 
quantity @ as (cf. flU), ©) 



- (TT^MWol*!) , 

(E-H + i7 1 )- 1 d 1 V . 



(50) 



Here E = E + e. Like the method of integral transforms, this procedure allows avoiding 
summation over final states at the calculation of a response function. However, in contrast to 
the state (0), or (|30|), the state (|50"D is not localized, and at large distances it includes outgoing 
waves in all open channels. The continuum-state calculation techniques are necessary to find 
this state which requires more effort then obtaining the localized state (p6|). The method was 



successfully applied for calculating three-nucleon responses [25.26 



VII. FEW NUCLEON RESPONSES FROM THE LORENTZ TRANSFORM 

These responses are smooth functions with single peaks and without pronounced additional 
structures which facilitates calculations. In Ref. [jTj] the longitudinal three-nucleon (e, e') re- 
sponse functions were studied. Eq. (|30] (for 0\ = O2) was cast into the form of inhomogenious 
Faddeev-like equations, and the source term was split into the three parts in a way preserving 
the permutational symmetries of Faddeeev components. Similar to Ref. 0, inversion was per- 
formed separately for the monopole contribution and for the sum of all other contributions to 
the response functions. The responses thus obtained compared well with those obtained from 
a conventional lengthy calculation, however the second above mentioned piece of the response 
function exhibited a small unphysical wiggle at low energy. When the wiggle was removed "by 
hand" the comparison with the conventional calculation became perfect. Probably one could 
have come to the same results if the proper threshold behavior of that piece of the response 
((e — e*/i) 3//2 in the 3 H case) had been imposed. It was demonstrated that the final results are 
nearly independent of 07, as they should. 

In Refs. [|§HT3"|| the three- and four-nucleon responses were obtained using an expansion 
technique in the dynamic equation ( |3Q]) {0\ = 2 ). The 4 He(e,e') longitudinal response func- 
tion M, the total 4 He 0, and tri-nucleon [TO-O, photodisintegration cross sections, and 
the 4 He spectral function along with the corresponding approximate 4 He(e, e') response [13 



were studied. Both central semi-realistic p-TO,0] and realistic j|TT|,|i"2"| NN interactions supple 



mented with a 3N force [12| were used. In [p|- [ll|JT3[] the correlated hyperspherical expansion f27| 



was applied to solve the problem. The basis functions were taken as products of the Jastrow 
factor, properly symmetrized hyperspherical harmonics coupled with the corresponding spin- 
isospin functions, and hyperradial functions. In Ref. [12|] the basis function are constructed as 



above mentioned products of spin-isospin-hyperspherical and hyperradial functions to which 
a Jastrow type spin-isospin dependent operator is applied. The prescription for selection of 
basis hyperspherical harmonics |[28|| , used in older work with no Jastrow correlations, is applied 
in conjunction with these correlations and is found to be very effective both for bound states 



and at solving Eq. (BO) [12 



It was found that the calculation of the transforms with an accuracy at a per cent level 
provides stable response functions. Achievment of such an accuracy in the transforms proved 
to be not requiring more effort than the calculation of binding energies with a comparable 
accuracy. As an example, in Fig. 2 the dependence of the Lorentz transform on the maximal 
value of the hyperspherical number K retained in the calculation is shown in the case of the 
four-nucleon photodisintegration response function. The calculation was done with the semi- 
realistic Malfliet-Tjon potential. The transforms at K max = 7 and at K max = 5 are already 
very close to each other. The sum rule checks (f44|) were satisfied by the responses obtained 
with a per cent accuracy. At the same time, let us mention that all possible checks of inversion 
results should be done, while the criterion of intermediate stability of the results with respect 
to the number N of basis functions in the expansion (^) may prove to be insufficient by 
itself [[11]]. Convergent 3-nucleon photodisintegration results with a realistic NN interaction 



supplemented with a 3N force are presented in [ ]T2"| . 

Finding practical ways to perform the calculations for larger systems, applying the approach 
to exclusive processes and to responses with more complicated structures, taking into account 
the long-range Coulomb interaction for exclusive reactions remains for future. 

The author wishes to thank W. Leidemann and G. Orlandini, his collaborators in the 
development of the present approach. The remarks on discretized responses were initiated by 
common calculations with I. Thompson. The work was supported by the Russian Foundation 
for Basic Research (grants no 96-15-96548 and 97-02-17003). 
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FIGURES 



The d(e, e') longitudinal response function. Solid curve: the conventional two-nucleon 
calculation. Dashed curve: calculation from the inversion of the Lorentz integral transform 
(07 = 5 MeV). 

The Lorentz transform L(a = —07? + 207(07?)) of the 4 He photodisintegration response func- 
tion at various K max values. The calculations were done for the Malfliet-Tjon NN potential. 
The 07 value ranged from 5 MeV to 20 MeV as a R increased. 
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